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to be related to previously observed phase transitions in the equilibrium 
flow topology. We find in these turbulent flow equilibria, two mechanisms 
for the appearance of ensemble equivalences, that were not observed in any 
physical systems before. These mechanisms are associated respectively with 
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mation) are discussed. It is found that decreasing the screening length scale 
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1 Introduction 

One of the major interests and achievements of equilibrium statistical me- 
chanics is the prediction of phase transitions: without describing details of 
the system configuration at a microscopic level, statistical mechanics predicts 
drastic changes in its macroscopic properties when varying only a few key 
parameters. Properties and computations of these phase transitions can be 
very different depending on the long-rangedness or short-rangedness nature 
of the interaction potential characterizing the physical system. We call long- 
range the interaction for which the potential is non-integrable. In the case 
of an algebraic potential V(r) ~ l/r°, this occurs when a is less than the 
(spatial) dimension of the system. 

Classical examples of long-range interacting systems are two-dimensionaQ 
and geophysical flows [55], self-gravitating stars and galaxies [TfllTSlIM] . or 
unscreened plasmas [H] . All these systems have the property of self-organizing 
spontaneously into large-scale coherent structures with a simultaneous com- 
plex small-scale dynamics: self-gravitating systems form galaxies and clusters, 
plasmas are organized into clouds of particles and turbulent flows create vor- 
tices and jets. Taking advantage of the long-rangedness of the interactions, 
statistical mechanics approach is particularly suitable to study these systems: 
it should predict the macroscopic organization of the system as the most 
probable outcome of its complicated microscopic dynamics, with constraints 
provided by the dynamical invariants, see e.g. the Miller-Robert-Sommeria 
equilibrium statistical theory (MRS theory hereafter) in the context of two- 
dimensional and geophysical turbulent flows [501EH1I5T1I52"] . 

Phase transitions and their related thermodynamic properties are well 
understood and have been thoroughly studied in the context of short-range 
interacting systems [33] ■ By contrast, despite the important progress in dy- 
namics and thermodynamics of long-range interacting systems achieved over 
the last decade (see e.g. [25JJ23 , 60 , 5 , 221I121IB"] ) and despite many analytical ex- 
amples of phase transitions in such systems (see e.g. [17] for two-dimensional 
flows, [2Z] for plasmas, [371I571US] for self-gravitating systems, [3D] for spin 
models), a clear understanding of phase transitions in the framework of long- 
range interacting systems is still lacking. 

One difficulty and challenging aspect associated with phase transitions 
in long-range interacting systems is the possible inequivalence between sta- 
tistical ensembles: a solution in the microcanonical (constrained) ensemble 
is not necessarily a solution in the canonical (unconstrained) ensemble [25j . 
The physical reason for this phenomenon is the non additivity of the energy 
E ~ J dr r d V(r) (with d the spatial dimension): when dividing a system into 
two subsystems, the energy of the whole system is not equal to the sum of the 
energy of each subsystem, because of the non-integrability of the interaction 
potential V . Then, many classical results of thermodynamics, obtained in the 
framework of short-range interacting systems, and relying on the additivity 
of the energy, are no more valid: i) statistical ensembles are not necessar- 
ily equivalent, ii) the equilibrium entropy S{E) may be non-concave. This 



The interaction potential is logarithmic for two-dimensional turbulent flows. 
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implies that the heat capacity C = (dT/dE) = -T- 2 d 2 S/dE 2 may be 
negative in the microcanonical ensemble: the temperature decreases when 
the energy increases. Such a counter-intuitive property was first predicted in 
the context of astrophysics by Lyden-Bell and by Thirring 59,3H133]- 

The appearance of ensemble inequivalence and negative heat capacity is 
associated with a zoology of phase transitions, with properties that would 
not be possible in systems with short-range interactions or local interactions 
only (for instance, first-order phase transitions with negative jump in tem- 
perature in the microcanonical ensemble). Such peculiarities have first been 
revealed in astrophysical context 00, 59 ,32 ,33.. Further characterization of 
ensemble inequivalence and its associated phase transitions have then been 
possible by analytical studies of toy models, as for instance (among others 
models) simple spin systems, see e.g. [21120] and references therein. Subse- 
quently, a classification of all possible phase transitions and their link with 
the occurrence (or not) of ensemble inequivalence when varying an external 
parameter was given [5]. Some of the transitions predicted have not been 
observed yet, neither in toy models, nor in real physical systems: this is the 
case for bicritical points (bifurcation from a first order to two second-order 
phase transitions) and for second-order azeotropy (simultaneous appearance 
of two second-order phase transitions) . 

Experimental observation of negative heat capacity and its related ther- 
modynamic peculiarities remains an open problem. Given a physical system, 
three necessary steps to tackle this problem are i) to predict a range of 
parameters associated with the occurrence of ensemble inequivalence, ii) to 
determine the nature of associated phase transitions iii) to relate these phase 
transitions to qualitative changes in the systems organization. Such is the 
focus of this paper, in the case of two-dimensional and turbulent geophysical 
flows. 

Remarkably, any piece of progress achieved for a specific long-range inter- 
acting system is useful for another one, since the computations of statistical 
equilibrium states generally involve similar variational problems. Let us for 
instance consider a physical system characterized by the mean field p at 
a macroscopic level, with one constraint given by the energy conservation 
£[p] = E. Assuming that the mean-field approach is valid, which is proven to 
be true for a large class of long-ranged interacting systems (see e.g [TrjIl2"5ll3l)I 

for the point-vortex model, [33] for two-dimensional Euler equations, and 
[25lfTj for a larger class of models), microcanonical and canonical problems 
can be described by two dual variational problems, respectively, 

S(E) = max.{S[p] | £[p] = E} and (1) 
p 

FQ3) = max [T[p] = -S[p] + p£\p]} , (2) 
p 

where S is a Boltzmann (mixing) entropy functional and T a (dual) Helmholtz 
free energy functional. These variational problems are dual because they 
have the same critical points, but not the same constraints. In presence of 
additional constraints, the only difference is the introduction of additional 
statistical ensembles associated with additional dual variational problems. 
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The validity of the mean-field approach is due the fact that the dynamics 
of the system at a particular location is not governed by local interactions, 
but rather by the whole system configuration: because every (infinitesimal) 
part of the system contributes to the local field, fluctuations around local 
mean-field values are small. 

The aim of the present work is to present analytical computations of 
phase diagrams for a large class of equilibrium states of two-dimensional and 
geophysical turbulent flows, in the framework of the MRS theory, emphasiz- 
ing the occurrence of ensemble inequivalence and its related thermodynamic 
properties. 

The interests are fourfold. 

First, beyond the application to two-dimensional and geophysical flows, 
the method used to compute the statistical equilibria is very general: in 
practice, we solve variational problems of the form ([1]) and ([2]), when the 
energy £ and the entropy S are quadratic functionals, and when an additional 
linear constraint is taken into accounlQ. Taking advantage of classical results 
about constrained and associated relaxed variational problems (a solution of 
([2]) is necessarily a solution of (jTJ , see for instance [25]), the computation of 
statistical equilibria will amount to demanding that a quadratic functional 
to be positive-definite. This computation will lead to an explicit criterion for 
the existence of statistical equilibria in different ensembles. 

Second, the computation of these equilibrium states relates previously 
observed phase transitions to the occurrence of ensemble inequivalence, thus 
deciphering its nature and physical origin. On the one hand, the existence of 
ensemble inequivalence in the context of two-dimensional flows was proven 
mathematically for point- vortices |11|I36| some time ago, numerically ob- 
served in quasi-geostrophic models [35] , as well as in a Monte Carlo study of 
point vortices in a disk [56]. On the other hand, phase transitions in the flow 
topology of equilibrium states have been described by Chavanis and Som- 
meria in the context of two-dimensional Euler flows [17] . The method used 
in the present paper to compute statistical equilibria allows for the simulta- 
neous description of these phase transitions and the appearance of ensemble 
inequivalence associated to them. In addition, we report in this context the 
first observation of a bicritical point and second-order azeotropy in a long- 
range interacting system (see also the recent works of [2T1I58] ). and discuss 
the related ensemble. 

Third, the class of models we study includes a 1.5-layer quasi-geostrophic 
model, which allows to discuss changes of phase diagrams in the presence 
of a screening length scale (Rossby radius of deformation). When this pa- 
rameter tends to zero, the energy is additive at leading order. However, it 
will be shown that non-local first-order corrections of the energy are essential 
to the computation of the equilibrium states, thereby explaining the appar- 
ent paradox that decreasing the screening length scale widens the ensemble 
inequivalence area in the phase diagram. 

2 In the context of two-dimensional flows, it will be shown that this corresponds 
to a class of solutions characterized by a linear relation between vorticity and 
streamfunction. 
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Fourth, a generalization of these results allows to describe a variety of 
"Fofonoff modes" characterized by a flow structure differing from the clas- 
sical low-energy solution of inertial academic ocean models [3D] ■ We provide 
a complete description of these modes, and explain why high-energy states 
could not have been obtained in previous computations of Fofonoff flows, for 
they did not take into account the possible existence of ensemble inequiva- 
lence. 

Some of these results have been already presented in a letter [61]. We 
present in this longer version detailed computations, discussions, and gener- 
alizations of this previous work. 

The paper is organized as follows: In section 2 we review a class of models 
for two-dimensional and geophysical turbulent flows, as well as the statistical 
mechanics of these systems. In section 3 we recall classical results on the 
equivalence between statistical ensembles, and explain how they can be used 
efficiently to compute equilibrium states of two-dimensional and geophysical 
flows; we then apply these methods to explicit calculations of the phase 
diagrams associated with any flow characterized by a linear relation between 
potential vorticity and stream function (q — ij) relation). In section 4 we 
present the thermodynamic properties of these systems first in the case of the 
Euler dynamics for two-dimensional turbulent flows; we report the generic 
existence of ensemble inequivalence, and describe phase transitions in the 
flow structure associated with it, including a bicritical point. In section 5 
we discuss generalization of these computations to a wider class of models 
and apply these results to a complete description of Fofonoff flows. This 
generalization makes it possible to describe second-order azeotropy. 



2 Equilibrium states of two-dimensional flows 

The aim of this section is to introduce a model for inviscid two-dimensional 
and geophysical turbulence (Euler equations and quasi-geostrophic equa- 
tions), and to review the equilibrium statistical mechanics of these systems. 



2.1 Models: Euler and quasi-geostrophic equations 

We consider in this paper a class of models that can be expressed as the 
transport of a scalar quantity q(r) by an incompressible, two-dimensional 
velocity field, where r = (x, y) is the spatial coordinate: 

d t q + u ■ Vq = with u=^-^, ^0 and q = A^~-^ + h, (3) 

where ip(j) is a streamfunction, h(r) a prescribed two-dimensional field in- 
terpreted as a topography (in addition, and without loss of generality, it will 
be assumed in the following that J v dr h = 0), and R the Rossby radius of 
deformation, an intrinsic length scale of the system. One must in addition 
specify a boundary condition for ip: we consider in this paper the case of 
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a closed domain T>, with impermeability constraint: there is no flow across 
the boundary dV. This means that the streamfunction is constant along the 
boundary. We consider here the simplest case ip = 0. Generalization to differ- 
ent ways of specifying the constant streamfunction at the boundarjU or even 
to other domain geometries (channel, doubly periodic) is straightforward. 

When R = +00 and h = 0, q is the vorticity and equations © are the two- 
dimensional Euler equations [22] • When R = +00 with possibly h 7^ 0, this is 
the barotropic quasi-geostrophic dynamics in presence of topography, which 
has been thoroughly studied in oceanic or atmospheric context |48) . and in 
which case q is the potential vorticity. When R is finite, this is the 1.5 layer 
quasi-geostrophic model (also called equivalent barotropic), a refinement of 
the barotropic model which includes the effect of stratification. This model 
is also relevant to describe plasma turbulence, in which case R is interpreted 
as the Larmor radius. When the topography is h = by, this is the Charney- 
Hazegawa-Mima model for either geophysical quasi-geostrophic turbulence 
on a beta plane or drift plasma turbulence [CTl31| . 

Conserved quantities. 

According to Noether's theorem, there is a dynamical invariant associated 
with each symmetry of the system, see e.g. |54j for a detailed discussion 
in the framework of geophysical flows. These dynamical invariants provide 
important constraints for the admissible states. 

The constraint associated with time invariance is the conservation of the 
total energy: 

m = \ ((w) 2 + ^) = -\ <(« -hw , (4) 

where (•) stands for a spatial integration over the whole domain T>, of area 
\T>\ = 1, and where the last equality is obtained by performing an integration 
by parts. 

Because the dynamics is the transport of potential vorticity q by an in- 
compressible flow, it conserves any Casimir functional 

CM = (.9(9)} , 

where g is any continuous function on T>. These conservation laws are equiv- 
alent to the conservation of the global distribution of potential vorticity, and 
are related to the particle relabelling symmetry [35]. In the present study, a 
Casimir functional that will be of particular interest (because of its simple, 
linear form), is the integral of the potential vorticity over the whole domain: 

C[q] = (q) . (5) 

3 In geophysical context, it is more physically relevant to consider ip = ipf r at 
boundaries, with ij)j r is determined by a mass conservation constraint reading 
f v dr xp = 0. The result presented in this paper would not be different in that 
case, but this lead to slightly more complicated computations. 
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If the domain has particular additional symmetries, for instance trans- 
lational invariance (case of a channel), or rotational invariance (case of a 
circular geometry), then additional conserved quantities have to be taken 
into account. We do not consider these situations in the present paper. 



Long-range interactions 

The energy (|4|) is the sum of a kinetic term (h ^(Vf/>) 2 ^ and a gravita- 
tional potential term ^ (^\\, but it can formally be expressed as a single 
potential energy term 

S[q] = -\ ((q h) ij) =-\J (q(r) - h(r)) G(r, r') (q(r>) - h(r')) dvdv' , 

(6) 



where Q (r, r') is the Green function solution of 



AQ 



g_ 

R 2 



with Q = for r on dT> 



(7) 



Far from boundaries, the green function is Q — j^Ko (l r — r '\/R) where Kq 
is the modified Bcsscl function of order zero. At short distance (|r — r'| <C 
i?), the potential vorticity reads q ~ Aip + h, and the dynamics is that 
of the 2D Euler equations (or of the barotropic model if h ^ 0), so the 
interaction potential has the classical logarithmic form: Kq ~ — In |r — v'\/R. 
At large distance (|r — r'| ^> R), the interaction potential becomes Kq ~ 

2 | r !L r /| J e~' r ~ r " R and is screened at scale R. If the screening parameter 

is much smaller than the typical length scale L = \/\T>\ of the domain where 
the flow take place (R <C L), then the system is short-ranged. Otherwise 
the system is long-range (for R = O(L) or R ^ L ). The screening length 
scale R is thus a parameter that allows one to pass from a short-range to a 
long-range interacting systems. 



2.2 Statistical mechanics of two-dimensional and geophysical turbulent flows 

The class of models we study are known to develop complex vorticity fila- 
ments at finer and finer scales. This makes almost impossible any attempt to 
have a deterministic approach of these systems, i.e. to describe the temporal 
evolution of the fine grained structure of the flow. Rather than describing 
these fine-grained structures, equilibrium statistical theories, assuming er- 
godicity, predict the final organization of the flow on a coarse grained level 
[521I531PH] . To keep track of the conservation laws of the transport equation, 
one needs to introduce the probability density p(r, a) of finding the potential 
vorticity level a at point r = (x,y), with a local normalization constraint 



N[p] = [ da pa = 1 , 
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where E stands for the admissible set of potential vorticity levels. The con- 
servation of the global vorticity distribution of potential vorticity is expressed 
as 



The energy of the microscopic potential vorticity fluctuations is negligible 
with respect to the mean- field energy of the statistical equilibrium state [52] , 
which can be rigorously justified using large deviation theory |42j . 

A macroscopic or coarse-grained state is then fully described by the prob- 
ability density p. The microscopic states are the fine-grained potential vor- 
ticity field of the flow. The number of microscopic states (any fine-grained 
field q) associated with a given macroscopic state (the probability density 
field p) can be quantified by the Boltzmann-Gibbs mixing entropy 



This has been rigorously justified using large deviations tools [151125) . Large 
deviation theory allows in addition to prove that an overwhelming number 
of microstates is associated with the most probable macroscopic state, i.e. 
the one that maximizes the mixing entropy while satisfying the constraints 
of the problem. To conclude, MRS statistical equilibria are solutions of the 
variational problem 



Smrs (E,d{&)) = max {S M rs [p] I £\p] = E & D a [p] = d(a) & N[p] = 1} , 



where E and d(a) are prescribed values of energy and given distribution of 
potential vorticity levels. Smrs (E, d (tr)) is the equilibrium entropy. Because 
this variational problem takes into account all the invariants of the dynamics, 
we call it full microcanonical variational problerr0. 

Computing the solutions of this variational problem (j5J) requires the 
knowledge of the infinity of Casimir constraints (or equivalently the knowl- 
edge of the distribution d(a)). This is a huge practical limitation to apply 
this theory, for two reasons. First, one does generally not know precisely what 
is the initial condition. Second, and more importantly, a variational problem 
with an infinite number of constraints is very difficult to handle. Many ef- 
forts have recently been devoted to alternative approaches which could lead 
to practical and mathematical simplification of this problem, see e.g. [231H], 
We present one of these methods in the following, and apply it to the actual 
computation of MRS equilibrium states. 

4 The term microcanonical will be used in the following when energy and circu- 
lation are the only constraints taken into account 





where tp is the solution of 






(8) 
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2.3 Simplification of the variational problem of the MRS theory 

Following the idea that less constrained, dual variational problems are easier 
to solve than more constrained ones (see e.g. [25]), it has been proposed to 
treat the Casimir invariants canonically to compute MRS equilibria [3], by 
studying the following variational problem: 

S(E, r) = max {S[q} = (s(q)) | £[q] — E k C[q] — -T} . (9) 

where s(q) is a concave function. Any solution to © is a MRS equilibrium, 
i.e. a solution of (|8]), but the converse may be wrong. A proof and a wider 
explanation of this statement is developed in 0]. The important point is 
that considering the coarse-grained potential vorticity fields q as an order 
parameter rather than the full distribution function (and treating canonically 
the constraints on the potential vorticity distribution) simplify the solution 
of the variational problem, not modifying the structure of critical points. 

We call this problem microcanonical by analogy with usual problems in 
thermodynamics, involving a constraint on the energy, and on the number of 
particle (here the circulation) in the maximization of a functional. Note that 
since a maximizer que of © is a coarse-grained potential vorticity field, say- 
ing that this solution is a MRS statistical equilibria means that there exists 
at least a set of constraints (the energy and the global vorticity distribution) 
for which the MRS equilibrium states Pmrs satisfy J da pmrs — QmC- 

The full microcanonical ensemble is relevant to describe an isolated sys- 
tem, in which the infinite number of constraints is fixed. In practice, this 
applies to situations where the inertial time scale of turbulent organization 
of the flow is much smaller than the dissipation and forcing time scales. 

We emphasize that considering a simpler variational problem, in a dual, 
less constrained statistical ensemble has no other reason that simplifying the 
mathematics. We interpret all the states that we compute as MRS equilib- 
rium state in the full microcanonical ensemble. 

A physical meaning of the other, less constrained ensembles would require 
the introduction of a bath of potential vorticity distribution, which would 
lead to fuzzy discussions, and has no reasonable physical justifications. 

Similarly, because the Casimir functional S[q] is maximized in the micro- 
canonical problem ©, it will be referred to as an entropy in the following, 
but only the Boltzmann-Gibbs (or mixing) entropy functional Smrs of the 
variational problem © has the physical meaning of an entropy. 

2.4 Critical points of the MRS variational problem 
2.4-1 Statistical equilibria are dynamical equilibria 

Critical points of the microcanonical variational problem satisfy SS — f35£ — 
jSr = 0, where the first variations are taken with respect to the potential 
vorticity field q, and where /J and 7 are the Lagrange multipliers associated 
with the energy and the circulation conservation respectively. This yields 
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g = s'- 1 (-^ + 7) , (10) 

where s /_1 is the inverse of the derivative of s, and where the inverse tem- 
perature /3 and the fugacity 7 are prescribed by the constraints on E and 
r. We conclude that critical points of (|9|) are flows characterized by a q — ip 
relation directly related to the function s. Such flows satisfy u • Wq = and 
are therefore dynamical equilibria of the transport equation ([3]) . 

2.4-2 The case of the maximization of a quadratic functional 

In the following, we shall compute solutions of the microcanonical problem 
© in the case of a quadratic functional by considering s(q) = — \q 2 : 

S[q\ = -\{?) , (11) 

The functional S[q] is the opposite of the enstrophy in that case. 

There are two reasons for studying this particular functional. First, it is a 
class of genuine statistical equilibria, with the property of being technically 
very simple. This example allows for a direct computation of the variational 
problem Q and gives analytic examples of ensemble inequivalence. Second, 
whatever the function s(q) is, computation of statistical equilibria in the limit 
of low-energy require studying the quadratic problem (fTTj) . 

Let us be more precise, considering the simple case when h = 0. The 
energy E is then a positive-definite functional of q. The limit of vanishing 
energy implies then that q is small. A multiplication of s(q) by a positive 
constant and the addition of a term linear in q does not change solutions of 
©. Assuming s twice differentiable, one can therefore consider without loss 
of generality that s(0) = 0, s'(0) = 0, and s"(0) = —1. A Taylor expansion 
gives then 

S (q) = -\q 2 + 0(q 3 ) . (12) 

In the case of a purely quadratic functional, using equation (|10[) . we find 
that the corresponding q — ip relation is linear: 

q = Pi* - 7 ■ 

This linear relation depends on two parameters, f3 and 7. These two param- 
eters have to be computed by considering the two independent constraints 
of the problem ©, namely the energy E and the circulation r. 

In the context of the Euler equations (R — +00, h = 0), Chavanis and 
Sommeria have computed and classified MRS equilibrium states correspond- 
ing to a linear q — ip relation, in an arbitrary closed domain |17) . They found 
two classes of phase diagrams, depending on the domain geometry. When 
the domain has a symmetry axis and is sufficiently stretched in the direction 
perpendicular to this axis, a transition from a monopole to a dipole is ob- 
served when varying the parameter E/T 2 above a critical threshold. When 
the domain is not stretched enough along the direction perpendicular to the 
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symmetry axis, no such transition is observed when varying the parameter 

E/r 2 . 

In [T7], the method to compute the MRS equilibrium states was, first 
to compute all the stationary solutions associated with a linear q — ip re- 
lation, and then to compute their entropy in order to select the entropy 
maxima. This made possible the observation (through analytical calcula- 
tions) of the intriguing phase transitions described above, but the precise 
nature of these transitions is still unknown, as well as the underlying phys- 
ical mechanism responsible for it. We show in the following that the use of 
another analytical method, taking advantage of the aforementioned relation 
between constrained and dual, less constrained variational problems, allows 
for the first time i) to compute explicitly statistical ensemble inequivalence 
in two-dimensional turbulent flows, ii) to relate this statistical ensemble in- 
equivalence to the occurrence of the phase transitions mentioned in previous 
paragraphs iii) to generalize these results to a larger class of models. 

3 Computation of MRS equilibrium states and ensemble 
inequivalence 

In this section we solve the variational problem © when S[q] = — (q 2 ) /2, 
and classify the solutions according to their energy and circulation. 

3.1 General method 

The strategy to solve © is to introduce and solve dual variational problems 
for which one or two of the constraints are relaxed. The natural first step 
is to consider the less constrained problem, referred to as grand- canonical, 
which is the easiest to solve: 

J(/3, 7 ) = mm{J[q] = -S[q] + /3 £[q] + 7 C [q]} , (13) 

where J is a thermodynamic potential. We compute its solutions, their energy 
and circulation. We then check which part of the plane E, r is filled by these 
solutions. Then we use the central result that any solution q gc of the relaxed 
variational problem ([13")) is also a solution of the more constrained variational 
problem ([9]) with constraints given by E = £[q c ] and r — C[q gc \. If some range 
of admissible parameters E, r are not achieved by grand-canonical solutions, 
one must consider a more constrained ensemble, and so on. In our case, it 
will be sufficient to consider the problem with only one constraint on the 
circulation, referred to as the canonical variational problem, in order to find 
all the solutions on the plane E, _T: 

Ftf, r) = min {T[q} = -S[q] + £[q] \ C[q] — T} , (14) 

9 

where F is a free energy. The fact that any value of the parameter (E, r) 
will be achieved by a canonical solution proves the equivalence between this 
ensemble and the microcanonical one. We will show in addition the existence 
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of a range of parameter E, r for which grand-canonical and canonical prob- 
lems are not equivalent: for this range of parameters, there is a canonical 
solutions, but no grand-canonical solutions. 

In the present case, the functional £ and S are quadratic and the func- 
tional C is linear. It is possible to take advantage of the linearity of this 
constraint in order to express the canonical problem (|14[) in the form of an 
unconstrained variational problem involving the minimization of a quadratic 
constraint. Then both problem (|13p and (|14l) are solved by finding the mini- 
mum of a quadratic functional, with a possible linear part and no constraints. 
Let Q and L be the linear operators associated to the purely quadratic part 
and to the purely linear part of this functional respectively. Then we have 
three cases: 

— Case 1. The smallest eigenvalue of Q is positive: the minimum exists and 
is achieved by an unique minimizer. 

— Case 2. At least one eigenvalue of Q is strictly negative. There is no 
minimum. 

— Case 3. The smallest eigenvalue of Q is zero (with eigenmode e neu t r ai)- 
If L[e neu trai] — (case 3a), the minimum exists, and each state of the 
neutral direction {ae neutra i\ QeK is a minimizer. If L[e neutra i} ^ ( case 
3b), then no minimum exists. 

We see that in the present case, the computation of statistical equilibria 
involves only the diagonalization of a quadratic operator. Application of this 
general method and detailed computations are provided in the remaining of 
this section. 



3.2 Decomposition on Laplacian eigenmodes 

In order to compute eigenvalues of the quadratic operators involved in the 
grand-canonical problem (|13[) and in the canonical problem (|14j) . it is con- 
venient to project the different fields (q, h and ip) on Laplacian eigenmodes. 
We introduce the complete, orthonormal basis {ei(r)} - gN of Laplacian eigen- 
modes on the domain T> 

(A - R- 2 ) e t = -(Me, = - (A, + R~ 2 ) * . (15) 

The Laplacian eigenvalues Xi are all positive, in increasing order. 

Potential vorticity q and topography h can be decomposed on this basis, 
as well as the streamfunction, computed by inverting the relation q — h = 

q = ^ Qjej h^^2h l e i if} = ^ — — —ei. 

Mi 

i i i 

The functionals S, £ and C given by equations (jllj) . and (J5]) respectively 
can then be expressed in terms of the coordinates {qi}'- 
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S[q) 



&4 



(16) 



£[q] 



(17) 



C[q] 



(18) 



where (e,) = f v dr e,(r) . 

We introduce two independent subspaces: 

— the subspace of the Laplacian eigenmodes having zero mean value ( (e£) = 
0, ieN*)and 

— the subspace of the Laplacian eigenmodes having non-zero mean value ( 



The notation "prime" and "double prime" will be used to distinguish (when 
necessary) one subspace from the other. 

The subspace of zero mean Laplacian eigenmode {e^} is generically empty 
when the domain geometry T> admits no particular symmetry: a small per- 
turbation of the domain would change the mean value of the eigenmode. If, 
by contrast, one imposes a symmetry axis in the domain geometry, then it 
exits generically eigenmodes having a zero mean: these eigenmodes are anti- 
symmetric with respect to the symmetry axis. This is for instance the case 
for a rectangular domain, for which the eigenmodes and eigenvalues are given 
in Annexe C. 

We assume that the smallest Laplacian eigenvalue on a closed domain is 
not degenerate, and that the corresponding eigenmode is positive everywhere, 
which turn out to be the case for simply connected bounded domains. It is 
thus an eigenmode of non-zero mean value: e\ = e", and /Uj = fi'{. Then, we 
have necessarily /i" < 

3.3 Solution for the grand-canonical problem 

We consider the grand-canonical problem (TT3"1) . Using expressions (fTB|) . (fTT|) 
and (|18p for entropy, energy and circulation functionals, we find: 



Since the quadratic part of J is diagonal in the Laplacian eigenmode ba- 
sis, we can immediately see that there is a unique solution to the variational 
problem if and only if j3 > — /ii (case 1. of subsection 13. ip . whatever the 
value of 7. If f3 < (case 2. of subsection 13. 1[) . there is no solution to 

the variational problem. If = —fii, then Q[ei] = (case 3. of subsection 



{e'O ± 0, % G W). 




(19) 
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13. ip and there is a neutral direction if and only if L[ei) = 0, which yields 
7 = hi I {ex). 

We show in Appendix [Blthat each point in the plane (E, r) located on the 
set of parabolas {Ep(r) \ /3 > —fix} corresponds to a unique grand-canonical 
solution, where Ef}(T) is given by 

w ,^ W + W1 r + (^g^)r*, (20) 

(see ([23])), and where ^[/i] = Bp[h] = for h = 0, and with 

/(/^Ettttt ( 21 ) 

i>X ^ l + p 

Interestingly, The curvature of the parabola Ep{r) equation (j2"0j) does 
not depend on the topography h, and is a decreasing function of (3. Notice 
also that the lower bound for the (admissible) energies among all flows with 
a given circulation is achieved when the inverse temperature tend to infinity 
j3 —> +oo. By taking this limit in equations (12T)|) . (l2~Tj) and using (h) = (see 
section 2.1), we find that these energy minima are located on the parabola 

E m (r) = B m [h]r + \(j2 ^ e *A r * > ( 22 ) 

where B m [h] = for h = 0. 

We conclude that grand-canonical solutions cover the whole area above 
the parabola E m (r) , given by (f2"2"j) and below the parabola E-^ (P) obtained 
by taking (i — > —fix in (|2T)]) (see also equation ([33]) in Appendix |B|) . 

Energies located above the parabola £L ^ (J 1 ) are not achieved by grand- 
canonical solutions, and the values located below E m (r) are not admissible. 
Because it exists a range of admissible energies and circulations not reached 
by grand-canonical solutions, we are in a situation of ensemble inequivalence. 
We then turn to the more constrained canonical problem to find solutions in 
this area. 



3.4 Solution for the canonical problem 

We now consider the canonical problem HU). It is possible to transform 
this constrained problem into an unconstrained variational problem, taking 
advantage of the linearity of the circulation constraint. Using this constraint, 
and recalling that (ei) ^ 0, one coordinate can be expressed in terms of the 
others: 

~ r — ^2 i>2 it i e i) 

q = q-qxex(r) with qx = j^t (23) 

\ex) 
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This expression is injected into the functional T = —S + (3£ appearing 
in the variational problem 



2 i>2 V f^iJ 7^2^ 



ij">2 



2 .4^. V y V /i'/y V Mi'/ (e\') 2 



1 1 



'I J I 



i>2 



e'l) A , /3 \ /3 t „ , (eft /8 



i / 



r±T^. 1 + 57 - 7^' + 77^5X' 8? , (24) 



where we have made the distinction made in subsection [32] between the two 
(independent) subspaces of Laplacian eigenmodes. The problem is now to 
find the minimizer q of this functional, with no constraint. 

To find the solutions to this problem, we use the same method as in the 
grand-canonical case. We call Q and L the linear operators associated with 
the purely quadratic and linear parts of T . 

In the subspace of zero-mean Laplacian eigenmodes, Q is diagonal. Its 
smallest eigenvalue is strictly positive if and only if j3 > — where \J X is 
the smallest eigenvalue of the linear operator — (A — Rr 2 ) associated with 
an eigenmode with zero mean-value. 

In the subspace of non zero-mean Laplacian eigenmodes, Q is not diago- 
nal, so this case requires more computations. We look for the value /3 = — /i* 
such that the smallest eigenvalue of Q is zero in this subspace (that corre- 
sponds to case 3. of subsection 13. We thus look for the eigenmode 

e* = e*ei such that Q\e*] = , 

i>2 

which yields: 

Let us first assume ^2j>2( e j) e j = ®- Equation (|25D then implies that 
there exists an integer k > 2 such that fj,* = fx^, with e* = for i ^ k, and 
e* k = a, where a a real. But then we have 2j>2( e i) e i = a ' which contradicts 
the initial hypothesis if a ^ 0. We conclude that ^2j >2 i e j) e j 0. 

Using this result, and multiplying (|2"5"j) by /Xi(ej)/ (m — fi*) yields then 

W \ 1 * A ^i( e i) 

vt > 2, e i = A , 

Mi - M 

where A can be determined by using a normalization condition for the eigen- 
mode. In addition, multiplying (|2"5|) by (ej)/ (1 — H* / Hi) and summing on 
i > 2 yields 

Mi(ei) 2 ^ Mi - M* 
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Finally, we obtain a condition that must satisfy ji* and the expression of 

e*: 



with e*^Aj2 J ^^^, (26) 



The main result is then that that — /i* is the smallest zero of the function^ 
/ given by equation (|21l) . Note that the class of solutions associated with 
j3 = —fj,* are q = ae* , with a any real number. Using equation (1231) . we find 
that the corresponding potential vorticity field are 

g = - — -ei + ae with e = e =- — ; ei (27) 

(ei) (ei) 

We conclude that whatever the value of r, there is a single solution 
to the variational problem (TH| if and only if (3 > — min{/i*, (case 1. 
of subsection 13. ip . When j3 < — min{/x*,^}, there is no solution to the 
variational problem (case 2. of subsection 13. X[) . When (3 = — min {/i* , \i\ } 
(case 3. of subsection 13. ip . we have to consider two cases depending on the 
sign of fi* — in order to discuss the existence of a neutral direction: 

— For /3 = — /z*, with /i* < h'i- Then a neutral direction with minimizers 
exists (case 3a of subsection 13. 1[) if L\e*] = 0. This leads to 



3^ V ( e i) 2 V~ Mi 7 Mi"' ' (ei) M 

By using (j2"6")l and a straightforward manipulation of the previous equa- 
tion, it yields the condition 

r*= M *£-^%. (28) 

For r ^ r* , there is no minimizer (case 3b of subsection 13. ip . 
— For (3 — —fi'i, with fi'i < fi* . Then a neutral direction with minimizers 
exists (case 3a of subsection 13.11) if L[e[] — 0. This gives the condition 
h[ = 0. If h[ ^ 0, there is no minimizer (case 3b of subsection 13. II) . 

From this analysis, we conclude that three different cases have to be consid- 
ered to describe the solutions in the phase diagram (T, E) , depending on the 
sign of fi'i — /i* and on the value of h[ (the projection of the topography on 
the smallest zero-mean Laplacian eigenmode) : : 



- i) n* < Hi or 

- ii) h* > Mi an d h[ = 

- iii) fi* > Mi an d h[ ^ 0. 



5 Using 1 = £j>i ( e i) e i and averaging this expression on the horizontal gives 
Z)i>i ( e *) 2 = 1 Z)i>i ( e *) 2 = !• Plugging tnis expression in fl2T} yields /(/3) = 
— 1 + /9 5^i>i ( e ») 2 / (^t + which corresponds to the function (3.8) obtained by 
Chavanis and Sommeria [17| . using another method. 
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For each of the three cases, the solutions of the canonical variational problem, 
as well as their energy and circulation, are computed in Appendix In all 
cases we find that all admissible values of circulation r and energy E are 
reached by these canonical solutions. We conclude the canonical problem is 
equivalent to the microcanonical one. 

In addition, since there is no grand-canonical solutions above the parabola 
E-^ t [T) given by equation (|35[) . grand-canonical and canonical ensembles are 
not equivalent for this range of parameters. 

The detailed computations of equilibriums states carried in Appendix [Q 
also lead to three different classes of phase diagrams, corresponding to the 
three cases above: 

— case i (/i* < there is a single equilibrium state at each point (T, E) 
of the phase diagram, except on the half line (r*, E > E*) where each 
point is associated with two states parameterized by two different values 
of 7. The parameters r* and E* are given by equations (|28I36[) . At high 
energy, the flow structure is dominated by the eigenmode e* defined by 
(US and (EH). 

— case ii > n[ and h[ = 0): there is a single equilibrium state at each 
point (r, E) below the parabola E_ fJ/ ' i (r) given by equation ([3"5]l . and two 
equilibrium states at each point (r, E) above this parabola. These two 
equilibrium states are parameterized by the value of their projection on . 
This value is zero on the parabola _B_ A1 ' i (/ n ), and tends to infinity with 
increasing energy. High-energy states are dominated by the eigenmode 
e[, which is a dipole with a neutral line along the symmetry axis of the 
domain. 

— case iii (//* > /Uj and h[ ^ 0): there is a single minimizer at each point 
(J 1 , E), and no peculiar canonical transition line. At high energy, the flow 
structure is dominated by the eigenmode e' x . 

The aim of the next sections is to detail an illustrate the peculiar thermo- 
dynamic properties of the system for these three classes of phase diagrams. 
In particular, we will clarify the nature of the phase transitions occurring on 
the half line F* , E > E* in case (i) phase diagrams, and on the parabola 
E_ fi > i (r) in case (ii) phase diagrams. Notice that these transition lines all 
stand in the ensemble inequivalence area, i.e. above the parabola .EL^ (r) 
described by equation ([35]). 



4 Phase diagrams in the Euler case 

For the sake of simplicity, we restrict ourself in this section to the case of the 
Euler equation (with h = 0, R — !► +00 in equation ([3])) in order to present 
the main striking features of the phase diagrams, namely the role played by 
the domain geometry, the changes in the flow structure associated with phase 
transitions, and the observation of a bicritical point. 
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Fig. 1 Phase diagram in case i (for a rectangular domain of aspect ratio 
r = 1.1 < r c ). Examples of flow structures are given on panel b. The ensemble 
inequivalence area is located above the blue parabola E\ 1 (r). The red line at 
r = is a first-order transition line, associated with a positive jump of 7. 



4.1 A geometry governed criterion 

In the case of the Euler equations, we have q = Atp, so the eigenvalues fa 
defined equation (fl~5|) are simply the Laplacian ones: fa = Aj. According to 
the previous section, there are only two types of phase diagrams when h = 0. 
The criterion for one phase diagram or another depends solely on the sign 
of fi* — Hi = A* — \' x , where A* is the smallest zero of the function (j2~Tj) 
Ei ^i( e i)/ C^t — A*) = 0) and where X[ is the smallest eigenvalue associated 
with a zero-mean Laplacian eigenmode. 

This criterion depends only on Laplacian eigenvalues, which depends 
themselves only on the domain geometry. We distinguish case (i) diagrams 
for which (A* < A^) from case (ii) diagrams for which ( A* > A^). 

If the domain admits no symmetry axis, only case (i) is possible, since 
there is generically no zero-mean Laplacian eigenmodes in that case. If the 
domain admits a symmetry axis, the sign of X[ — A* must be computed. In 
the case of a rectangular domain, there is a critical aspect ratio that can be 
computed numerically. This computation gives t c ~ 1.12, where the aspect 
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Fig. 2 Phase diagram in case ii (for a rectangular domain of aspect ratio 
r = 1.8 > r c ). Examples of flow structures are given on panel b. The ensemble 
inequivalence area is located above the blue parabola E\ 1 (r). The Green dashed 
parabola Ey (F) is a second-order transition line. Inside this parabola, d 2 S/dr 2 > 




ratio is defined by r = L x /L y with L x and L y the lengths of the rectangular 
domain), and which was already reported by Chavanis and Sommeria j!7j . 
For aspect ratios smaller than r c , the phase diagram is in case (i), since 
Aj — A* > 0. For aspect ratios greater than r c , the phase diagram is included 
in case (ii), since \[ — A* < 0. It is expected (but not proven) that any 
domain geometry admitting a symmetry axis, and sufficiently stretched in a 
direction perpendicular to this axis, is in case (ii). 

The figures presented in this section are computed in the case of a rect- 
angular domain, in which case the Laplacian eigenvalues and eigenmodes are 
explicitly known (see Appendix [A~|) . According to the analysis of the previous 
section, these results are however generic to any domain geometry. Since all 
the two-dimensional fields (<7, h, ip) are decomposed on Laplacian eigenval- 
ues in the present study, the explicit computation of these quantities (and 
of the associated functionals £[q], C[q] and S[q]) in an arbitrary geometry 
would require in practice a truncation of this decomposition, and the use of 
a numerical solver to compute Laplacian eigenmodes and eigenvalues in the 
given domain geometry. 
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4.2 Thermodynamic properties of the phase diagrams 

From the knowledge of the Laplacian eigenvalues, one is able to draw the 
equilibrium entropy S(E, r) in terms of the internal parameters E and r, 
using the computations carried out in the previous section (the expression of 
the entropy in terms of the projections qi is given by (jlip . and the expressions 
of these projections g, are given by equations (|2T)1) and (1571) ). 

We present such plot of S(E, r) on figures [l}c (for case (i) phase dia- 
grams) and[2]-c (for case (ii) phase diagrams). All the thermodynamic prop- 
erties of the phase diagrams can be deduced from these plots. 



Ensemble inequivalence area. 



The most striking feature appearing on figures [T]-c and ®c is the existence 
of a region of parameters E, r for which S(E, r) and its concave hull do 
not coincide. Following the analysis of the previous section, this region of 
parameters is the one above the parabola 

£_ Al (r) = — 4-t r 2 , 

2Ai(ei) 

obtained by taking h = and \i\ = Ai in equation (|35[) . This parabola is 
represented as a blue line on figures [TJa and [2}a, and corresponds to the 
upper boundary for grand-canonical solutions. This illustrates the result of 
[25j : there is ensemble inequivalence whenever the entropy does not coincide 
with its concave hull. 

This is to our knowledge the first analytical description of such a generic 
occurrence of ensemble inequivalence in two-dimensional flows. Chavanis and 
Sommeria observed (using analytical computations) phase transitions in the 
flow structure for the Euler equations, but without relating them to the 
appearance of ensemble inequivalence |17j . Ellis, Haven and Turkington did 
compute numerically a specific case of statistical equilibria showing ensemble 
inequivalence, but without understanding of what control the appearance 
of this situation, and without deciphering the link between the ensemble 
inequivalence and the phase transitions |26] . 

We describe in the following peculiar microcanonical phase transitions 
associated with the existence of such ensemble inequivalence. 



Phase transitions for case (i) phase diagrams. 



By taking h = in equations ([25 ]) and ([36 ]) . we obtain r* — 0, E* = 0. The 
half line 

(r = Q,E > 0) 

is represented as a thick red line in figure [T}a. We observe on panels [T]-c 
and [TJ-e the occurrence of a discontinuity of 7 = dSJ dT across this line (see 
Appendix iDl for explicit computations). This is the signature of a first-order 
phase transition. We also note that this discontinuity corresponds to a posi- 
tive jump of 7. Such a thermodynamic peculiarity would not be possible in 
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an ensemble equivalence area: a positive jump of the derivative of the en- 
tropy (dS/dr) can only occur in a region where the entropy S(E, E) and its 
concave envelope do not coincide. 

Whatever the values of the energy E and of the circulation r, there is a 
single equilibrium state, having the structure of a monopole, as represented 
at low- and high-energies on figure [TJd, except on the first-order transition 
line, represented in red on figure QJa. When the first-order transition line 
is crossed, the flow structure changes from a monopole of a given sign to 
a similar monopole with the opposite sign, with coexistence of both states 
on the transition line. Notice that the high-energy states tend to a solution 
dominated by eigenmode e* (when f3 — » —A*), given by equations (|26p and 
(|27p . The structure of this eigenmode can be computed numerically. In the 
case of a rectangular domain of aspect ratio r < 1.12, it is always a monopole. 

Phase transitions for case (ii) phase diagrams. 

We have shown in the previous section that above the parabola E_y (E) 
given by equation with V% fa = A*, 

^(/(A^f^-Ai) 2 ; 

there are two entropy maxima (depending of the sign of their projection on 
the eigenmode e^) for each point E, E of the phase diagram, while it exists 
only one such entropy maximum for each point of the phase diagram below 
this parabola. This parabola is represented as a green line on figure [U-b. 

When the parabola E_x[ (E) is crossed, there is a discontinuity of d 2 S/ dE 2 
and of d 2 S/dE 2 , as observed on figures [2}c,e and as explicitly computed in 
Appendix [Dj This is the signature of a second-order phase transition. 

Notice the existence of positive values of d 2 S/dE 2 in figure [5J This last 
property is similar to negative heat capacity, but for the parameter E rather 
than for the energy E. This corresponds to a region of parameters E, E for 
which the entropy S(E, E) is convex: it does not then coincide with its con- 
cave envelope. We conclude that such a peculiarity can occur only in the 
ensemble inequivalence area. 

Below the parabola E_y {E), there is a unique equilibrium state associ- 
ated with each point (E, E). These low-energy states have the structure of a 
monopole. 

Above the parabola -E-a; (E), each point corresponds to two equilibrium 
states. For a given circulation, the different states above the parabola differ 
only by the value of their projection on e\ (a dipole). The choice of one 
state among the two possibilities above the parabola breaks the system's 
symmetry. At high energy, this contribution dominates: the flow is therefore 
a dipole. 

For a fixed circulation, there is thus a continuous transition from a monopole 
(at low-energy) to a dipole (at high energy). 



22 




r 



Fig. 3 Bicritical point: bifurcation from a first-order phase transition line (in 
red) to two second-order phase transitions lines (in green). Case of a rectangular 
domain, for the 2-D Euler equation. The aspect ratio r is taken as an external 
parameter, the circulation r as an internal parameter, and the energy is fixed. Small 
insets are schematic representations of the entropy S [q] in which the two directions 
correspond to projection onto the eigenmodes e* (horizontal) and e'i (vertical). M 
, m and s stand respectively for M aximum, m inimum and s addle 



Observation of a bicritical point 

Let us fix the energy and vary the aspect ratio r in the case of a rectangular 
domain. Our previous analysis predicts a transition from case (i) to case (ii) 
phase diagrams, above a critical value t c of the aspect ratio. The correspond- 
ing phase diagram (r, r) is presented in figure |3] There is a bifurcation from 
a first-order transition line (corresponding to the red line in case (i) phase 
diagram plotted in figure [1} to two second-order transition lines (correspond- 
ing to the green plain-dashed double line in case (ii) phase diagram plotted in 
figure[5]), at the point t c = 1.12, r = 0. Such a bifurcation from a first-order 
to two second-order transitions is referred to as a bicritical point. 

Bicritical points have already been observed (and actually defined) in the 
context of short-range interacting systems, see e.g. |29j . However, there is to 
our knowledge no example of such bifurcations for systems with long-range 
interactions, while their possible existence was predicted in [5]. More impor- 
tantly, the previous analysis shows that this bicritical point occurs in the 
ensemble inequivalence area, and is in that respect a signature of ensemble 
inequivalence, associated with drastic changes in the flow structure (a tran- 
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Fig. 4 Schematic representation of ensemble inequivalence, for a fixed energy and 
a functional 5 = (s(q)) with s(q) = s(—q). Inserts represent S[q] for a given value 
of the parameter F. there are two symmetric solutions for _T = 0. 



sition from a monopole to a dipole). 

Small insets on figure [3] are schematic representations of the functional 
S[q] = -<q 2 )/2forafixed energy and different values of the parameters r, r. 
The horizontal direction corresponds to the eigenmode e* defined by (|2l)|) and 
(l27l) and associated with eigenvalue A* . The vertical direction corresponds to 
the eigenmode e^, the smallest Laplacian eigenmode with zero mean- value, 
associated with A^. 

Below the green dashed line (corresponding to ft — —X^) the maxima of 
S[q] have no contribution along e[, consistently with the computation of the 
corresponding equilibrium states, given by (1291) . 

On the the line r = 0, for aspect ratio r < r c , there are two entropy 
maxima proportional to e*, with j3 = —A*. (When varying circulation, one 
of these two maxima becomes metastable.) 

On the line r, for aspect ratio r > r c , there are two entropy maxima 
proportional to e' 1; with ft = —X'i- 

The bicritical point corresponds to the transition from case (i) to case 
case (ii), for which A^ = A*. Equilibrium states are degenerate at this point, 
since any linear combination of the two corresponding eigenmodes (satisfying 
the constraint) are entropy maxima. 
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4.3 Fixed energy and varying circulation: generic occurrence of ensemble 
inequivalence, analogy with a magnetic system 

Strikingly, in all the phase diagrams described above, the ensemble inequiv- 
alence, visualized when the entropy does not coincide with its concave hull, 
is visible when drawing S(E, E) with E fixed and r varied, but not visible 
when drawing S(E, r) with r fixed and E varied. This phenomenon were 
already mentioned in 26 . Similarly, we report positive values of d 2 S/dE 2 , 
but not of d 2 S/dE 2 (that would correspond in that case to negative heat 
capacity). 

One can understand why ensemble inequivalence occurs generically for 
the parameter r in variational problems of the form © when C = (q), and 
with the symmetry s(q) — s(—q). Let us consider to simplify the case of a 
domain without symmetry, for a given, fixed energy. 

Then the symmetry q — > —q is associated with the generic occurrence 
of a first-order transition line for E = 0, on the plane EE. The reason 
for the occurrence of this phase transition is that they are necessarily two 
symmetric equilibrium states on the line E — 0: if q is a solution of the 
variational problem j9j, then — q is also a solution. These two states are on 
two branches of solutions, symmetric with respect to r = 0, see figure 2] 
Varying the circulation E at fixed energy E favors one of the configurations 
(the other becoming metastable) , which leads to a microcanonical first-order 
phase transition. 

Using then the fact that first-order phase transition in the constrained 
(microcanonical) ensemble generically implies inequivalence with the less con- 
strained (canonical) ensemble (see figure 2 of [5]), we conclude that there 
is necessarily ensemble inequivalence for the parameter E around the line 



This case of ensemble inequivalence is associated with the breaking of 
the discrete symmetry q —q. It would similarly occur in a magnetic sys- 
tem with long-range interactions, for which the computation of statistical 
equilibria would lead to solve the problem: 



where m is the magnetization, h is an imposed magnetic field, and where 
e(m) = e(— m) and s(m) = s(— m). Considering a fixed energy etot such 
that the equilibrium states have a non zero magnetization when h = 0, 
and considering the magnetic field h as an external parameter, there is an 
ensemble inequivalence area associated with the occurrence of a first order 
microcanonical phase transition at h = 0. In that respect the circulation E 
of a two-dimensional flow plays the role of an external magnetic field h in a 
magnetic system. It has to be noted however that this symmetry breaking is 
partial as the "field" E acts only on the eigenmodes with < e, >^ 0. 



r = o. 



S{e to t, h) = max \ s(m) | e tot = e(m) + hm 



m 
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5 Generalization and application to Fofonoff flows 

We explain in this section how the phase diagrams of the Euler case are 
changed when a non zero topography or a finite Rossby radius are taken into 
account. Considering the effect of the Rossby radius of deformation allows for 
a discussion of the behavior of the system in presence of a scale that screens 
the interactions. Considering the effect of topography can breaks some of 
the system symmetries, and allows for the description of a new peculiar 
thermodynamic properties, namely second-order azeotropy. It also allows for 
a description of Fofonoff flows in the framework of RMS statistical theory. 



5.1 Effect of the topography h(x,y): symmetry breaking and second-order 
azeotropy 

We consider the effect of a topographic term h in equation ([3]) by considering: 

q = Aip + h , 

which has two physical consequences: i) a symmetry breaking, ii) the possible 
presence of second-order azeotropy. 

5.1.1 First physical consequence: symmetry breaking 

First, a term h(x, y) generically breaks the symmetry q — > —q. Consequently, 
phase diagrams in parameter space (P, E) will generically be non-symmetric 
with respect to the axis r = 0, unless one imposes a peculiar symmetry to 
h. 

Second, even if the domain geometry admits some symmetry, the topo- 
graphic term h generically breaks this symmetry, unless one imposes this 
symmetry to h. This symmetry breaking is associated with the possibility of 
a third class of phase diagrams, referred to as case (iii), when the domain has 
a symmetry axis and is sufficiently stretched in the direction perpendicular to 
its axis. We have seen previously that in case (ii) phase diagrams, i.e. when h 
satisfies the symmetry of the domain geometry, there is a second-order phase 
transition line in the parameter space (r, E), above which the system admits 
two possible equilibrium states, differing only by the sign associated with the 
mode e[ (a dipole). In case (iii) phase diagrams, the topographic term favors 
one of these two degenerate equilibrium states (i.e. it selects the sign of the 
dipole e'i), and there is therefore no second-order transition line. 

To conclude, in presence of topography, there are three classes of phase 
diagrams. The criteria for these phase diagram are explicitly computed in 
appendix [Cj as well as the corresponding thermodynamic properties. The 
result depends only on the value of the projection {he'i) of the topography 
on the eigenmode e[ : either h[ = (case with symmetry) , or h\ =^ (generic 
case without symmetry), and on the the eigenvalues m defined equation (fT"5)) . 
which, in the present case (q = Atp + h), are simply Laplacian eigenvalues 
(/i. t = Xi): 
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i) (A* > A'jJ. There is a first-order micro-canonical transition half line, as in 
the Euler case, but the minimum of this line (the point (J 1 *, E*) defined by 
equations ([25)) and (I5rj)) ). is located strictly above the parabola E^\ 1 (r) 
given by equation (|55|) . 

ii) (Ai > A* and h[ — 0). There is a second-order transition line as in the 
Euler case, namely the parabola E_x[(r), given by equation The 
curvature radius of the parabola does not depend on h, but the minimum 
of the parabola does, and it is now located strictly above -E_Ai(-0- This 
case corresponds to figure [5l 

iii) (X[ > A* and h[ ^ 0). There is no microcanonical phase transition above 
i?_Ai(^), because the topography breaks the symmetry satisfied by the 
domain symmetry. 

Note that in these three cases, there is still ensemble inequivalence above 
the parabola E-^^r) given by equation (|35l) . For h ^ 0, the minimum of 
this parabola is no more located on (_T = 0, E — 0), but its curvature radius 
is left unchanged. 

Note also that except its role in breaking the symmetry of the system, 
the topography does not induce qualitative changes in the phase diagrams 
described in the Euler case. As for the structure of the equilibrium states, the 
topography plays an important role for low-energy states, but high-energy 
states are those of the Euler dynamics. 

The energy is a quadratic positive-definite functional of the field (q — h). 
Low energy states are therefore those for which q ~ h: the structure of 
the potential vorticity field (and consequently of the streamline) is entirely 
determined by the topography field. 

By contrast, the contribution of the topography becomes negligible for 
high-energy states (for which q = Alp + h « Atp), which are therefore those 
of the Euler dynamics. 

5.1.2 Second physical consequence: second-order azeotropy 

A case (ii) phase diagram associated with a topographic term h(y) is repre- 
sented in figure[5]-a . Below the parabola E^\ 1 (r) represented as a blue plain 
line, the entropy is a concave function of the circulation, so 7 — dS/dT is a 
strictly decreasing function. When the minimum C of the parabola E_\ ± (r) 
is crossed with increasing energy, there is the occurrence of an ensemble in- 
equivalence area, associated with a first-order transition in the statistical 
ensemble with fixed energy and relaxed circulation. This first-order transi- 
tion is illustrated with the Maxwell construction figure[S]-c. We conclude that 
point C is a critical canonical point. 

When the minimum A of the parabola E_\ t ' (r) represented as a plain- 
dashed green line is crossed with increasing energy, there is the simultaneous 
appearance of two second-order phase transitions in the microcanonical en- 
semble. The discontinuity of dj/dT associated with these two second-order 
phase transitions can be seen in figure EJ-d, point F. 
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The simultaneous appearance of two phase transitions was referred to as 
second-order azeotropy in the classification [5]- The term azeotropic point is 
commonly used for binary mixture to define the point of the phase diagram 
where there is the simultaneous appearance of two first-order phase transi- 
tions. This term were originally used because there is no change in the frac- 
tional composition of a binary mixture when the azeotropic point is crossed. 
The term "second-order azeotropy" extends this definition to the outbreak 
of two second-order phases transitions from nothing. However, contrary to 
the case of first-order phase transitions in binary mixtures, there is no phase 
separation associated with these second-order phase transitions, but rather 
a superpositions of one phase on the other when the transitions occur: we 
have seen previously that the two phases involved in the second-order tran- 
sitions of case (ii) phase diagrams represent two different flow structures: 
a low-energy state superimposed with a "dipole phase" whose contribution 
increase with the energy. 

To our knowledge, second-order azeotropy has never been observed in any 
physical system with long-range interactions, excepted in the recent study of 
[58] in the context of the HFM model. First order azeotropy has also been 
reported for the first time only recently |21j . As explained in [S], there are 
two types of second-order azeotropy: one is visible in the canonical ensemble, 
the other is not. The azeotropic point we describe fall in this second class. All 
the thermodynamic properties related to the second-order azeotropy appear 
as predicted in [5]: 

— Sufficiently close to the azeotropic point A, the discontinuity of 7 is not 
associated with a change of concavity for the entropy: the slope of 7(I 7 ) 
does not change its sign. 

— For energies sufficiently larger than the energy of the azeotropic point A, 
the discontinuity of d^/dr is associated with a change of concavity of 
the entropy: the slope of j(r) changes sign, as for instance seen in figure 
[SJ-e, point H. 



5.2 Application to an academic inertial ocean model, the Fofonoff flow 

We now apply the results obtained previously to the description of an aca- 
demic ocean model known as the Fofonoff solution, also referred to as Fofonoff 
modes or flows. Fofonoff flows are stationary solutions of the undissipated un- 
forced barotropic quasi-geostrophic model in a closed domain £>, described 
by equation ([3]), with h = —by, R = +00 and q = (3tp — 7: 

Aip — fiip = by — 7 

with ip = on dT> . 

Fofonoff computed and described the solutions of these equations inde- 
pendently of statistical theories, assuming (3 ^> 1 [30]. The only reason to 
consider such an approximation was to be able to solve the mathematics. 
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Fig. 5 Second-order azeotropy: simultaneous appearance of two second-order 
phase transition (at point A). Case of a rectangular domain of aspect ratio r = 
L x /L y > 1.12 with h = sm(jiy / L y ). 



Within this approximation, the contribution of the Laplacian term is negli- 
gible in the domain bulk: 



It corresponds to a weak westward current (the velocity in the west-east 
direction is u x = —d y ^p ~ —b//3) taking place in the whole domain, and 
with strong boundary jets that close the circulation. The emergence of such 
state in freely evolving numerical simulations has been previously observed 
and discussed by [63 , 64 , 24 , and interpreted as statistical equilibrium states 
predicted by energy-enstrophy statistical theory in presence of topography, 
which leads to linear q — ip relations [55115] . 

Energy enstrophy theory has been first proposed by Kraichnan in the 
framework of the statistical mechanics of truncated Euler equation [38] . A 
drastic consequence of the truncation is that only the energy and the en- 
strophy are conserved quantities of the dynamics, while any higher vortic- 
ity moments are conserved for the Euler dynamics. This energy-enstrophy 
approach has then been generalized to a class of geophysical flows that in- 
clude the effect of bottom topography and of the variations of the Coriolis 
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Fig. 6 Fofonoff flows, for the three possible phase diagrams (i), (ii) 
and (iii) The circulation is fixed, slightly greater than zero, and h = by. The 
streamfunction is represented in a rectangular domain with different aspect ratios 
(cold and hot colors are negative and positive values respectively). At low energy, 
the flow is always the usual Fofonoff mode, with a weak westward flow in the 
domain bulk, and strong recirculating jets at the boundaries. If the domain is not 
stretched enough (case (i)), the high-energy state is a monopole. If the domain is 
sufficiently stretched in the East- West direction (case (ii)), there is a second-order 
phase transition: the high-energy state is one dipole in the East-West direction. 
The choice of sign of this dipole breaks the system symmetry. If the domain is 
sufficiently stretched in the North-South direction (case (iii)), there is no phase 
transition. The high-energy state is a dipole displaying the Fofonoff mode structure, 
but without strong jets confined at the boundaries. 



parameter with latitude, as well as the stable stratification of these flows 
[55] . This energy-enstrophy approach has been shown to be equivalent to the 
phenomenological minimum enstrophy principle of [5], that were proposed 
independently in the same geophysical context |13) . It has subsequently been 
shown that equilibrium states of these statistical mechanics approaches cor- 
respond to a particular class of MRS statistical equilibria [I], namely the 
class of equilibrium states computed in the present work. 

Only low-energy states of barotropic ocean models on a beta plane (such 
as the original solution of Fofonoff) were previously computed in the frame- 
work of statistical theories. Our results make possible a direct computation of 
all MRS equilibrium states associated with a linear q — ip relationship model, 
whatever their energy and circulation. The interest is twofolds: 

— First, it shows that the classical Fofonoff flows are MRS equilibria located 
in the ensemble equivalence area: there are low-energy states character- 
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ized by an inverse temperature j3 1, for which there exists correspond- 
ing to grand-canonical solutions. 
— Second, we compute also high-energy states, and find that such states can 
have a flow structure very different from the classical Fofonoff solution, 
when they are located in the ensemble inequivalence area, see figure [6] 

In addition, all the phenomenology related to phase transitions and ensemble 
inequivalence applies for those models: according to the domain geometry, one 
can obtain the phase diagram corresponding to case (i), (ii) or (hi) described 
in the previous sections (see figure [6]) . 

In particular, any domain that admits a symmetry axis in the y direction, 
and that is sufficiently stretched perpendicularly to this axis, corresponds to 
the phase diagram presented on figure [5l for which there is second-order 
azeotropy. 

5.3 The effect of the Rossby radius of deformation 

We now discuss the effect of a finite value for the Rossby radius of deformation 
R by considering: 



see equation ([3]) . The only change induced by this term is that the eigenvalues 
\ii (solutions of [A — i? -2 ) ej = — /Xi6j) are now Laplacian eigenvalues Xi 
shifted by — I / R 2 . The existence and structure of the phase diagrams (i) , (ii) 
and (iii) are left unchanged, since the criterion fi[ — fi* = X[ — X* does not 
depend on R. 

However, since the first Laplacian eigenvalue is of order L~ 2 , where L is 
a typical length scale of the domain, taking the limit R <C L in (|35| gives 
(r) ~ r 2 / (2R- 2 ) (ei) 2 . When R < L, the curvature radius of this 
parabola tends to infinity. Remembering that this parabola is the boundary 
for ensemble inequivalence, the range of parameters E, r associated with 
the inequivalence between canonical and grand-canonical ensembles tends 
therefore to fill the whole parameter space E, r. 

This result might at first sight seem surprising: the screening length scale 
R tends to zero, which means that the energy is an additive quantity at lead- 
ing order, and yet one observes an ensemble inequivalence area larger than 
in the Euler case ! We have explained that such inequivalence between sta- 
tistical ensembles can occur for long-range interacting systems which 
means that non-local terms must remain important in some way to set the 
structure of the equilibrium state, even if the energy is additive at leading 
order. 

Let us be more precise: writing q = £ + xJj/R 2 , with £ = Aip, the energy 
reads 



6 Ensemble ine quiv alence can also occur in finite volume, small short-interacting 
systems, see e.g. |18| . but not when the thermodynamics limit is considered. 
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In the limit R <^ L, the term (ip 2 ) / R 2 dominates over the term (£1/)), which 
is the one involving long-range interaction. Similarly, the enstrophy reads 




and is dominated by the term (ip 2 ) /R 4 when R <C L. Then, at leading 

order, S = —R~ 2 £, which suggests to write j3 = —R~ 2 + f3. Finally, the 
dominant term of £ and S cancel each other when computing —S + f3£. To 
find canonical solutions, one looks for the minimum of 

- S + ? £ =^( m + ^ + °{^)) ' 

with the constraint on the circulation. At leading order, this problem is then 
to find the minimum of (£ip) + (V' 2 ) with the constraint on the circulation; 
it is independent of R, and involves the nonlocal term {£ip}. This shows that 
this non-local term, even if negligible with respect to the total energy, is 
crucial to compute the statistical equilibrium state. 

The fact that the screening length scale plays no peculiar role in the 
structure of the statistical equilibria can also by seen by considering critical 
points of the variational problems, which satisfy q = /3 (ip — (ip))+r. Writing 
then j3 — f3 + R~ 2 gives 

py>-(rl>)) + r = Aii> , 

which is just the equation satisfied by statistical equilibria of the Euler dy- 
namics in the case of a linear q — tp relationship. An important physical 
consequence is also that the parameter R plays no particular role in the flow 
structure in that case. 

These results and peculiarities rely on the strong assumption of a purely 
linear q — ip relation (related to a purely quadratic functional S). Previous 
analytical studies of statistical equilibria in the limit R <C L, and in the 
case of tanh-like q — tp relationship, lead to very different flow structure 
and thermodynamic properties than in the linear case, with, for instance, 
no ensemble inequivalence, and the formation of strong jets of width R in 
the flow structure [3] . Such non-linear q — ip relation are the relevant one for 
geophysical applications, as for instance the explanation of Jovian vortices 
[5] or oceanic rings and jets [52] . 

More detailed computations and studies will therefore be necessary to 
fully understand the role played by the parameter R when non-linearities in 
the q — Tp relation are taken into account. This will be the focus of a future 
work. 



5.4 Generalization to a larger class of models 

The previous computations apply beyond the particular case of two-dimensional 
and geophysical flows. It can be extended to the computation of statistical 
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Fig. 7 Schematic representation of the functional — S[q] + f3£[q] in configura- 
tion space, where /3 is such that a) — S[q] + f3£ is negative-definite (only stable 
directions); b) — S[q] + fi£ has a neutral direction; c) — S[q] + f3£ has an unsta- 
ble direction. Only case b leads to non trivial solutions (the points of the neutral 
direction are represented as a dashed red line). 



equilibrium states of any system for which the energy is a quadratic func- 
tional, in the limit E — > 0. This could for instance be applied to the com- 
putation of statistical equilibria of the Vlasov equations. More generally, the 
results presented in the present paper hold true whenever one has to solve a 
variational problem that consists in maximizing a negative-definite quadratic 
functional I[q] = —S[q] of a field q(r), with constraints given by another def- 
inite positive quadratic functional £ [q] = E and a linear functional C[q] = r. 

If in addition the constraints are respectively purely quadratic (i.e. with 
no linear nor constant part) and purely linear (i.e. with no constant part), 
then there is generical ensemble inequivalence. For the sake of simplicity, let 
us (for instance) consider the line r — in the phase diagram (r, E): on this 
line, there are no grand-canonical solutions (i.e. no solution in the ensemble 
without constraints), but there are canonical solutions (i.e. solutions in the 
ensemble with only the circulation constraint. A schematic representation of 
the computation of such canonical solutions is given figure [JJ It sums up the 
general method followed in this paper to compute statistical equilibria: For 
j3 above some critical value, the (purely quadratic) functional — S + j3£ is 
negative-definite, and its minimum is located at q = 0. For j3 below some 
critical value, there is at least one unstable direction (one eigenmode with 
negative eigenvalue). For f3 equal to the critical value, there is a neutral 
direction, corresponding then to non trivial statistical equilibria: these are 
the solutions of the line r = in the phase diagram (r, E) . 

The important conclusion is the generic occurrence of ensemble inequiv- 
alence for this class of systems. Properties of phase transitions inside the 
ensemble inequivalence area depend on the properties of a system at hand 
(in particular the existence of particular symmetries) . 



6 Conclusion 

The main result of this paper is the observation of a generic ensemble in- 
equivalence region for a wide class of models, including two-dimensional and 
quasi-geostrophic turbulent flows, through explicit analytical computations. 
We have shown that the ensemble inequivalence is related to the occurrence 
of peculiar phase transitions, associated with drastic changes in the flow 
topology. Several possible phase transitions have been identified in parame- 
ter space (r, E), leading to three different kinds of phase diagrams. We have 
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shown that the criterion for these different phase diagrams depends on the 
domain geometry only, and can be explicitly computed for a system at hand. 

Strikingly, statistical ensemble inequivalence is generically observed for 
varying circulation at fixed energy; we have explained that this is related to 
existence of a symmetry go — q in the system when it admits no particular 
spatial symmetry. 

When varying the energy at fixed circulation, we have shown the exis- 
tence of marginal cases of ensemble inequivalence at sufficiently high-energy 
(entropy then becomes a linear function of energy). Taking into account non- 
linearity is necessary to discuss these marginal cases. We conjecture that 
tanh-like relations lead to ensemble equivalence (concave entropies), while 
sinh-like relations could be associated with ensemble inequivalence areas. 
Understanding the role of this nonlinearity will be the focus of a future work 

m 

Some of the transitions and related ensemble inequivalence we described 
were predicted in the classification [5], but they had not been observed in 
any physical system yet. This is the case of 

— Second-order azeotropy: the outbreak from nothing of two second-order 
phase transitions, associated with a switch from a monopole to a dipole, 
breaking the system's symmetry when the energy is increased and the 
circulation r is kept fixed. 

— Bicritical point: the bifurcation from a first-order transition line to two 
second-order transition lines when increasing the domain aspect ratio, 
at fixed energy, and considering the circulation r as the only internal 
parameter. The first-order phase transition is associated with a switch 
from a monopole to another one; the two second-order transition lines 
correspond respectively to a switch from a monopole to a dipole, and 
to a switch from this dipole to another monopole, symmetric (equal and 
opposite) to the initial one, when varying the circulation. 

Finally, we have pointed out a surprising effect induced by a screening 
length scale (the Rossby radius of deformation) in 1.5 quasi- geostrophic mod- 
els. Decreasing this length scale widens the ensemble inequivalence area in 
the phase diagrams while the interaction becomes more local. We have ex- 
plained this apparent paradox, showing that non-local first-order corrections 
of the energy are essential to compute equilibrium states, even if the energy 
is additive at leading order. We have shown that this phenomenon is peculiar 
to the class of equilibrium states considered in this paper (states character- 
ized by a linear q — ip relation). For such states, the flow structure does not 
depend on this screening length scale, which can be absorbed in the q — tjj 
relation. We conjecture that such peculiarity would not be observed in phase 
diagrams obtained by considering a varying energy and fixed distributions of 
potential vorticity levels. 

All states computed in this paper could have been observed in the frame- 
work of energy-enstrophy theories, which all lead to the same variational 
problem as the one we have studied in this paper. Most of the solutions 
to this variational problem were however not described in previous studies, 
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carried out in the framework of energy-enstrophy theories, because the com- 
putations were only performed in the grand-canonical ensemble, where all 
the constraints are relaxed. 

Kraichnan noted the existence of a " condensed states" for particular val- 
ues of the temperature, which correspond to the microcanonical states that 
we describe in the case r = (no circulation) [35]. He concluded that such 
condensed states were the most interesting ones, but that there was in av- 
erage no mean flow: for all values of the parameter E, there were always 
two equiprobable symmetrical states, cancelling each other in average. Con- 
sequently to these studies, it has often been said that energy-enstrophy the- 
ories do not predict the emergence of a mean flow when no topography is 
taken into account. The results of the present paper clarify this issue: when 
performing computations in the microcanonical ensemble, we have seen that 
there exist (possibly multiple) equilibrium states for each value of the circu- 
lation and the energy, so that a mean flow is possible. 

In the presence of topography, just as in the case of Fofonoff flows, com- 
putations in the relaxed ensemble lead to low-energy statistical equilibria, 
even with zero circulation. But just as in the Euler case, high-energy states 
could not be obtained in the most relaxed, grand-canonical ensemble. It is 
necessary to work in the canonical ensemble to describe them. This explains 
why only low-energy, Fofonoff modes, had been described by previous com- 
putations of statistical equilibrium states in the presence of topography, see 
e.g. [551]. 

Some of statistical equilibria were already computed and described in the 
Euler case by [17] , in the framework of the MRS statistical theory, but using 
another method, and without reference to thermodynamic properties of the 
phase diagrams in relating to the physics of long-range interacting systems. 
Our contribution, following |61j . has been to link these previously observed 
phase transitions to the generic occurrence of statistical ensemble inequiva- 
lence, using a method that can be applied to a class of systems much larger 
than two-dimensional and geophysical turbulent flows. Subsequently to our 
work, some of these results on the thermodynamic properties of phase dia- 
grams have also been presented in [46] (case without topography) and [45] 
(case with topography), using the method and formalism of |17j . 

One of the current challenging problems in physics of long-range inter- 
acting systems is to come up with experimental setups for which ensemble 
inequivalence properties such as negative heat capacity could be observed. 
So far, no physical system has been proposed to realize such an experiment. 
The results presented in this paper indicate that two-dimensional flows are 
serious candidates to first observe ensemble inequivalence. 

Phase transitions associated with the existence of ensemble inequivalence 
are predicted in the microcanonical ensemble. This ensemble can show up 
if the time scale for inertial organization of the flow is much smaller than 
the forcing and dissipation time scales. When such separation of time scale 
does not exist, it is more relevant to include effects of forcing and dissipation. 
The understanding of phase transitions taking place at equilibrium is then a 
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starting point for studies out of equilibrium [7] . 



A Laplacian eigenvalues of a rectangular domain 

We give here the Laplacian eigenmodes (Aei = — Xiei) in a rectangular domain. 
Instead of indexing the eigenmodes of the Laplacian by i, we use two indices m, n 
corresponding to the wave numbers in the x- and y directions, respectively: 

e mn — 2 sin (imx/yfr) sin (^my^fr) , 

where r = L x /L y and the domain is supposed to be of area unity (L y L x = 1). The 
corresponding eigenvalues are 

2 f n 2 2 \ 

Amn = 7T I — + Tin I , 

and {emn} = 8/ (nmiv 2 ) for n and m odd, (e mn ) = for n or m even. The zero- 
mean eigenmode associated with the greatest eigenvalue — Ai is ei = en and the 
zero-mean eigenmode associated with the grea test eigenvalue — X[ is e[ = ei 2 (if 
r > 1). We see when considering equation l|22C that in this case, E m (r) is the line 
E = in the plane (_B, r). 



B Computation of energies and circulations of grand-canonical 
solutions 

Let us compute the energy and the circulation of all the grand-canonical solutions, 
by considering separately the cases B > —fJ,i and /3 = —fii. In any case, a direct 
computation of the critical points of (I19|l gives 

Vi > 1, (fK + P)qi = -(ym(ei)+phi). 

For p > and an arbitrary value of 7, the critical points are unique and 
given by 

V»>1, *G9,7) = -^-t; (iM^i) + Phi) ■ (29) 
A straightforward computation gives their circulation: 

r = £9i<ei> = -7/(/?)-/3£-^rl (30) 
i i>i ^ + p 

where 

/(/ ? ) = E^t5- (31) 

Their energy E = \ J2i {qi — hi) 2 / fXi can be expressed either as a function of 
the Lagrange multiplier 7 

* w " (g itf&r) +1 (g (g • <»» 
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or as a function of the circulation r, by using expression (|30fl : 

Ea(T) = A a [h] + Ba[h]r + ( — y M*(e*) 2 \ p2 

where A a [h] = Ba[h] — when h = 0, whatever the value of f3. 

For ft = — Hi, the solutions of the neutral direction ( 7 = — /ii/(ei}) are param- 
eterized by a £ R: 

qi = a and Vt > 1, ©(0,7) = ( ^rhi + hA . (34) 

M» - Mi \< e i) / 

By computing the circulation _T of these states, we find that a — _T/(ei): for a 
given circulation, there is a single solution in the neutral directio n. I t is straight- 
forward to check that this solution corresponds to the solution (|30fl in the limit 
f3 — > —Mi- Then, the solutions of the neutral direction are located on the parabola 

e_ w (r) = A- n [h] + b_ w [h]r + 1 r 2 , (35) 

2Mi(ei)' ! 

obtained either by using the expression (|34[l or by taking the limit /? — > — mi in (|20l l. 

We conclude that each point of the phase diagram [E, r) located on the set 
of parabolas {Ea(r) j f} > — Mi) corresponds to a unique grand-canonical solution, 
where Ea{F) is given by (|33|) . 



C Computation of energies and circulations of canonical solutions 

The aim of this Appen dix is to compute energy and circulation of the solutions of 
the canonical problem (|14[1 . 

When P > — min {m*, m'i}i using equation (|30f) . one can express 7 as a function 
of r since /(/3) 7^ 0, and substitute this expression into equation (|29[) to obtain 
the exp ression of the solutions q(P,T). The corresponding energy Ea(r) is given 
by PD|. 

There is a single canonical solution at each point of the region spanned by 
the set of parabolas {Ea(r) | f3 > — min{Mi, M*}}- Since min{/ii, M*} > Mi) an d 
since the curvature of these parabolas is a decreasing function of /?, we can see that 
the set of grand-canonical solutions is included in the set of canonical solutions. 

Three cases must be considered to describe the solutions in the phase diagram 

(r,E), 

— i) M* < Mi or 

— ii) m* > Mi an d h'i = 

— iii) m* > Mi an d h'i 7^ 0. 



Case i ( fi* < 

There is a n eutr al direction for j3 — — m* if and only if _T = _T* (where F* is given 
by equation l|28|)). In that case, 7 can not be determined by equation (|30p . Without 
specification of the energy, 7 can then be chosen arbitrarily on R. By considering 
the expression (|29p we see that the parameter 7 parameterizes the equilibrium 
states of the neutral direction: 

37 (~M , " ) =7e H — , 
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where e* is given by equ ations (|26[) and (|27l) . The energy (7 ) of these states 

is given by equation (|32[) . It varies between a minimum value _E* and +00, where 

E*= min ( 7 )} . (36) 

Notice that for each point _T*, i? > there are two states parameterized 
by two different values of 7, and that if _T 7^ _T*, then the energy (|20|l diverges 
when ft — > — /i*. The set of parabolas {^(.T) | /3 > — ^i*} hence covers the half 
plane E, F above E m (F), except the half line (T = T*, E > E*) which is filled by 
equilibrium states of the neutral direction. 

Case ii (fi'i < fi* and h[ = 0). 

There is a neutral direction whatever the value of r, when f3 = —f/i. The states of 
this neutral direction are 

q = q(-fi' 1 ,r) + ae' 1 , a£R , (37) 
where q(—fi'i,r) is given by combining (|29|) and (|30p . Their energy is 

E = E_, S r) + ±- a \ 

where E_ fl i (r) is obtained by taking f3 = —fi'i in (|20[) : 

e., (r) = A-,, [h] + b_„, [h]r + [ - ) nn2 y - ^\ ) r 2 . (38) 

The set of parabolas {£^(F) | /? > — spans the whole range of parameters 
(E,r) located above E m {r) and below E^(r). The canonical solutions of the 

neutral direction, parameterized by F and a — q[ (the projection on e'i), are 

located on the set of parabolas j-E^ (F) + a 2 /(2fj,' 1 ) \ a £ r|, which covers the 

half plane above E^i (F). 

Case Hi ( ji'-y < fj,* and h[ ^ 0). 

In that case, there is no neutral direction. One can easily show that the energy of 
the minimizers obtained in the case /3 > —Hi diverges when j3 — > The set of 

parabolas \Ep(F) \ f} > —fi'i} covers the whole half plane above E m (r). 



D Computation of phase transitions in the Euler case. 

D.l First-order phase transitions. 

Let us compute the discontinuity of the entropy across the half line (F — 0, E > 0). 
We have 

s = -\( q 2 ) = -^{ q ^) + lr = pE + lr . (39) 

Let us consider the limit ft — > —A* for a given value of the energy. The relation 
(|30|l with n* = A* gives 

r = if (-A*) (0 + A*) (40) 



:;<s 



and then S = -X*E + -yr/2 + FEj (7/' (-A*)) + O ({P + A*) 2 ) . Similarly, the 
equilibrium state (|29p reads (with first-order correction) q — 'ye* + O (/3 + A*) , 
where e* is given by ((26]) and ([27]), with Vi, fn = Ai. Then E = -^7 2 (eM _1 e*) + 
O (/3 + A*). By noticing that (e*A~ 1 e*) = f (-A*) and by using 00]), we obtain 

7 = ±K*\/s + o (r/\/S) with K * = \J f,^} > ( 41 ) 

which gives g = ±K*\/E~e* + O \ T/\/~E \ . To determine the sign of the maximizer 

of S at fixed energy and circulation, we inject ()41|) in (|39p an use (|40|) to find the 
equilibrium entropy 

S = -A*E + Kl-/E|r| +0 

Thus, the sign of the maximizer q is determined by the sign of the circulation. 
When the circulation zero, both states ±g are maximizers. Since k* does not vary 
with r, there is a positive jump of 7 = dS/dF on _T = 0. 

We conclude that the half line r*,E > E* is a first-order transition line. 
Two states (parameterized by two different values of 7 = dS/dF) coexists at each 
point of this line, and any perturbation of the circulation select one of those states 
as the unique equilibrium state. The entropy is a concave function of the energy 
(d 2 S/dE 2 < for f / 0), with a marginal situation on the first-order transition 
line (d 2 S/dE 2 = for T = 0). 




D.2 Second-order phase transitions. 

Let us compute this second-order phase transition, by considering the limit f3 — > 
— Ai for a given energy. We have 5 = fiE + yF/2. Two cases have to be considered 
: a) P = —X'i b) < —X'i- In case a) (in which E,T are located above E_ X > 1 ), 
7 = r/f(—Xi) and S = — X^E + r 2 /f(—X' 1 ) whatever the circulation and the 
energy. Then d 2 S/dE 2 = and d 2 S/dr 2 = d'y/dr = 0. In case b), the relation 
Ptij) becomes 

The equilibrium state can be expanded at lowest order in /3 + X[: 
q = 7ed + 7 (P + Ai) e c2 + O (r (P + X[)) . 

with 

E-° Ai(ei) Aj(ej) 
— — ei and e c2 =^- — ^, 

»=i Al A i i= i l A * _ A J 

The energy reads £ = 7 2 £ c i+7 2 (/? + Ai) E c2 +0 {P + Ai) with E cl = -§ (ed^-^d), 
_E C 2 = — \ (e C 2A~ 1 e c i). By using the previous expression of 7, we obtain 

E = TTT^Twi ((1 - 2/' (- A '0 (0 + A '0) ^ + (P + A '0 ^) + {P + A '0 - 
which gives /3 as a function of the energy and the circulation: 



S c2 - 2£ c i/' (-Ai) 
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Then 

d 2 s _ dp _ ff(-K) 



8E 2 BE \ T J E c2 -2E cl f'(-X' 1 ) ' 

and 

which is generically not equal to zero when r 0. We see that d 2 S/dE 2 and 
d 2 S/dF 2 are discontinuous when the parabola E_y i (r) is crossed, i.e. when /3 = 

Since the second-order derivative of the entropy toward the toward the energy 
and the circulation are discontinuous when this parabola is crossed, we conclude 
that this is a second-order transition line. 



Acknowledgements We thank J. Barre, P.H. Chavanis, and J. Sommeria, Ana 
Carolina Ribeiro-Teixeira for interesting discussions, Marianne Corvellec for provid- 
ing detailed and helpful comments on the manuscript, and an anonymous reviewer 
for pointing to us the analogy with magnetic systems, for their comments on this 
manuscript. This work was supported by the ANR program STATFLOW ( ANR- 
06-JCJC-0037-01 ) and the ANR program STATOCEAN (ANR-09-SYSC-014). 



References 

1. Barre, J.., Bouchet, F., Dauxois, T., Ruffo, S.: Large Deviation Techniques Ap- 
plied to Systems with Long-Range Interactions. Journal of Statistical Physics 
119, 677-713 (2005). DOI 10.1007/sl0955-005-3768-8 

2. Barre, J., Mukamel, D., Ruffo, S.: Inequivalence of ensembles in a system with 
long-range interactions. Phys. Rev. Lett. 87(3), 030,601 (2001) 

3. Bouchet, F.: Mecanique statistique des ecoulements geophysiques. PHD, Uni- 
versite Joseph Fourier- Grenoble (2001) 

4. Bouchet, F.: Simpler variational problems for statistical equilibria of the 2d 
euler equation and other systems with long range interactions. Physica D 
Nonlinear Phenomena 237, 1976-1981 (2008) 

5. Bouchet, F., Barre, J.: Classification of Phase Transitions and Ensemble In- 
equivalence, in Systems with Long Range Interactions. Journal of Statistical 
Physics 118, 1073-1105 (2005). DOI 10.1007/sl0955-004- 2059-0 

6. Bouchet, F., Gupta, S., Mukamel, D.: Thermodynamics and dynamics of sys- 
tems with long-range interactions. Physica A Statistical Mechanics and its 
Applications 389, 4389-4405 (2010) 

7. Bouchet, F., Simonnct, E.: Random Changes of Flow Topology in Two- 
Dimensional and Geophysical Turbulence. Physical Review Letters 102(9), 
094,504 (2009) 

8. Bouchet, F., Sommeria, J.: Emergence of intense jets and Jupiter's Great Red 
Spot as maximum-entropy structures. Journal of Fluid Mechanics 464, 165- 
207 (2002) 

9. Bretherton, F.P., Haidvogel, D.B.: Two-dimensional turbulence above topog- 
raphy. Journal of Fluid Mechanics 78, 129-154 (1976) 

10. Caglioti, E., Lions, P.L., Marchioro, C, Pulvirenti, M.: A special class of sta- 
tionary flows for two-dimensional Euler equations: A statistical mechanics de- 
scription. Commun. Math. Phys. 143, 501-525 (1992) 

11. Caglioti, E., Lions, P.L., Marchioro, C, Pulvirenti, M.: A special class of sta- 
tionary flows for two-dimensional Euler equations: A statistical mechanics de- 
scription. Part II. Commun. Math. Phys. 174, 229-260 (1995) 

12. Campa, A., Dauxois, T., Ruffo, S.: Statistical mechanics and dynamics of solv- 
able models with long-range interactions. Physics Reports 480, 57-159 (2009) 

13. Carnevale, G.F., Frederiksen, J.S.: Nonlinear stability and statistical mechanics 
of flow over topography. Journal of Fluid Mechanics 175, 157-181 (1987) 



40 



14. Charney, J.G.: On a physical basis for numerical prediction of large-scale mo- 
tions in the atmosphere. Journal of Meteorology 6(6), 372-385 (1949) 

15. Chavanis, P.H.: Phase transitions in self-gravitating systems: Self-gravitating 
fermions and hard-sphere models. Phys. Rev. E 65(5), 056,123 (2002) 

16. Chavanis, P.H.: Phase Transitions in Self-Gravitating Systems. International 
Journal of Modern Physics B 20, 3113-3198 (2006) 

17. Chavanis, P.H., Sommeria, J.: Classification of self-organized vortices in two- 
dimensional turbulence: the case of a bounded domain. J. Fluid Mech. 314, 
267-297 (1996) 

18. Chomaz, P., Gulminclli, F.: Phase Transition in Small System. Nuclear Physics 
A 749, 3-13 (2005) 

19. Corvellec, M., Bouchet, F.: A complete theory of low-energy phase diagrams 
for two-dimensional turbulence equilibria. In preparation (2011) 

20. Costeniuc, M., Ellis, R.S., Touchette, H.: Complete Analysis of Phase Tran- 
sitions and Ensemble Equivalence for the Curie- Weiss-Potts Model. J. Math. 
Phys. 46, 063,301 (2005) 

21. Dauxois, T., de Buyl, P., Lori, L., Ruffo, S.: Models with short- and long-range 
interactions: the phase diagram and the reentrant phase. Journal of Statistical 
Mechanics: Theory and Experiment 6, 15 (2010) 

22. Dauxois, T., Ruffo Cugliandolo, L. (eds.): Long-range interacting systems: 
lecture notes of the les houches summer school: Vol 90, august 2008 (2009) 

23. Dauxois, T., Ruffo, S., Arimondo, E., Wilkens, M. (eds.): Dynamics and Ther- 
modynamics of Systems With Long Range Interactions (2002) 

24. Dukowicz, J.K., Greatbatch, R.J.: Evolution of Mean-Flow Fofonoff Gyres in 
Barotropic Quasigeostrophic Turbulence. Journal of Physical Oceanography 
29, 1832-1852 (1999). DOI 10.1175/1520-0485(1999)029 

25. Ellis, R.S., Haven, K., Turkington, B.: Large Deviation Principles and Com- 
plete Equivalence and Nonequivalence Results for Pure and Mixed Ensembles. 
J. Stat. Phys. 101, 999 (2000) 

26. Ellis, R.S., Haven, K., Turkington, B.: Nonequivalent statistical equilibrium 
ensembles and refined stability theorems for most probable flows . Nonlinearity 
15, 239-255 (2002) 

27. Euler, L.: General principles of the motion of fluids. Physica D: Nonlinear 
Phenomena 237(14-17), 1825 - 1839 (2008). Euler Equations: 250 Years On - 
Proceedings of an international conference 

28. Eyink, G.L., Spohn, H.: Negative-temperature states and large-scale, long-lived 
vortices in two-dimensional turbulence. Journal of Statistical Physics 70, 833- 
886 (1993) 

29. Fisher, M.E., Nelson, D.R.: Spin Flop, Supersolids, and Bicritical and Tetra- 
critical Points. Physical Review Letters 32, 1350-1353 (1974) 

30. Fofonoff, N.P.: Steady flow in a frictionless homogeneous ocean. J. Mar. Res. 
13, 254-262 (1954) 

31. Hasegawa, A., Mima, K.: Pseudo-three-dimensional turbulence in magnetized 
nonuniform plasma. Physics of Fluids 21, 87-92 (1978) 

32. Hertel, P., Thirring, W.: Free energy of gravitating fermions. Commun. Math. 
Phys. 24, 22-36 (1971) 

33. Hertel P. Thirring, W.: Soluble model for a system with negative specific heat. 
Annals Phys. 63, 520-533 (1971) 

34. Ispolatov, I., Cohen, E.G.D.: Phase transitions in systems with l/r a attractive 
interactions. Phys. Rev. E 64(5), 056,103 (2001) 

35. Kiessling, M.K.H.: Statistical mechanics of classical particles with logarithmic 
interactions. Comm. Pure Appl. Math. 47, 27-56 (1993) 

36. Kiessling, M.K.H., Lebowitz, N.: The Micro-Canonical Point Vortex Ensemble: 
Beyond Equivalence. Letters in Mathematical Physics 42(1), 43-56 (1997) 

37. Kiessling, M.K.H., Neukirch, T.: Negative specific heat of a magnetically self- 
confined plasma torus. Proc. Nat. Acad. Sci. 100, 1510-1514 (2003) 

38. Kraichnan, R.H., Montgomery, D.: Two-dimensional turbulence. Reports on 
Progress in Physics 43, 547-619 (1980) 

39. Landau, L.D., Lifshitz, E.M., Reichl, L.E.: Statistical Physics, Part 1 (Third 
Edition). Physics Today 34, 74-+ (1981). DOI 10.1063/1.2889978 



41 



40. Lynden-Bell, D., Wood, R.: The gravo-thermal catastrophe in isothermal 
spheres and the onset of red-giant structure for stellar systems. Mon. Not. 
R. Astron. Soc. 138, 495 (1968) 

41. Macor, A., Doveil, F., Elskens, Y.: Electron Climbing a "Devil's Staircase" in 
Wave-Particle Interaction. Phys. Rev. Lett. 95(26), 264,102-+ (2005) 

42. Michel, J., Robert, R.: Large deviations for young measures and statistical 
mechanics of infinite dimensional dynamical systems with conservation law. 
Communications in Mathematical Physics 159, 195-215 (1994) 

43. Michel, J., Robert, R.: Statistical mechanical theory of the great red spot of 
jupiter. Journal of Statistical Physics 77, 645-666 (1994) 

44. Miller, J.: Statistical mechanics of euler equations in two dimensions. Phys. 
Rev. Lett. 65(17), 2137-2140 (1990). DOI 10.1103/PhysRevLett.65.2137 

45. Naso, A., Chavanis, P.H., Dubrulle, B.: Statistical mechanics of Fofonoff flows 
in an oceanic basin. ArXiv e-prints (2009) 

46. Naso, A., Chavanis, P.H., Dubrulle, B.: Statistical mechanics of two- 
dimensional Euler flows and minimum enstrophy states. ArXiv e-prints (2009) 

47. Padmanabhan, T.: Statistical mechanics of gravitating systems. Phys. Rep. 
188, 285 (1990) 

48. Pedlosky, J.: Geophysical fluid dynamics. New York and Berlin, Springer- 
Verlag, 1982. 636 p. (1982) 

49. Ripa, P.: Symmetries and conservation laws for internal gravity waves. In: 
American Institute of Physics Conference Series, American Institute of Physics 
Conference Series, vol. 76, pp. 281-306 (1981) 

50. Robert, R.: Etats d'equilibre statistique pour Pecoulement bidimensionnel d'un 
fluide parfait. C. R. Acad. Sci. 1, 575-578 (1990) 

51. Robert, R.: A maximum-entropy principle for two-dimensional perfect fluid 
dynamics. J. Stat. Phys. 65, 531-553 (1991) 

52. Robert, R., Sommeria, J.: Statistical equilibrium states for two-dimensional 
flows. J. Fluid Mech. 229, 291-310 (1991) 

53. Robert, R., Sommeria, J.: Relaxation towards a statistical equilibrium state in 
two-dimensional perfect fluid dynamics. Phys. Rev. Lett. 69(19), 2776-2779 
(1992) 

54. Salmon, R.: Lectures on Geophysical Fluid Dynamics. Oxford University Press 
(1998) 

55. Salmon, R., Holloway, G., Hendershott, M.C.: The equilibrium statistical me- 
chanics of simple quasi-gcostrophic models. Journal of Fluid Mechanics 75, 
691-703 (1976) 

56. Smith, R.A., O'Neil, T.M.: Nonaxisymmetric thermal equilibria of a cylindri- 
cally bounded guiding-center plasma or discrete vortex system. Physics of 
Fluids B 2, 2961-2975 (1990) 

57. Stahl, B., Kiessling, M.K.H., Schindler, K.: Phase transitions in gravitating 
systems and the formation of condensed objects. Planet. Space Sci. 43, 271- 
282 (1995) 

58. Staniscia, F., Chavanis, P.H., De Ninno, G.: Out-of-equilibrium phase transi- 
tions in the HMF model: a closer look. ArXiv e-prints (2010) 

59. Thirring, W.: Systems with negative specific heat. Zeitschrift fur Physik 235, 
339-352 (1970). DOI 10.1007/BF01403177 

60. Touchette, H., Ellis, R.S., Turkington, B.: An introduction to the thermody- 
namic and macrostate levels of nonequivalent ensembles. Physica A Statistical 
Mechanics and its Applications 340, 138-146 (2004) 

61. Venaille, A., Bouchet, F.: Statistical Ensemble Inequivalence and Bicritical 
Points for Two-Dimensional Flows and Geophysical Flows. Physical Review 
Letters 102(10), 104,501 (2009) 

62. Venaille, A., Bouchet, F.: Solvable phase diagrams and ensemble inequivalence 
for two-dimensional and geophysical turbulent flows. ArXiv e-prints, submitted 
to Journal of Physical Oceanography (2010) 

63. Wang, J., Vallis, G.K.: Emergence of Fofonoff states in inviscid and viscous 
ocean circulation models. Journal of marine research pp. 83-127 (1994) 

64. Zou, J., Holloway, G.: Entropy maximization tendency in topographic turbu- 
lence. Journal of Fluid Mechanics 263, 361-374 (1994) 



